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One  of  the  most  important  methods  of  obtaining  laser  action,  both 
from  a  technological  and  scientific  viewpoint,  is  the  chemically- 
pumped  laser.    The  HN     +  CO     chemical  laser  has  been  modeled  and 
an  experimental  investigation  conducted  to  determine  the  existence  of 
laser  action.    The  analytical  model  has  been  used  to  predict    the  expected 
population  inversion  and  laser  gain  for  a  range  of  pressures  and  various 
additives.    The  experimental  investigation  conducted  for  comparison  with 
the  model  did  not  result  in  laser  action.    The  increase  in  pressure  for 
the  mixtures  investigated  indicates  the  occurrence  of  an  explosion.    The 
lack  of  detection  of  laser  action  has  been  attributed  to  difficulties  with 
alignment  of  the  optical  system  and  to  signal  noise  associated  with  the 
flash  circuit. 
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I.     INTRODUCTION 

The  first  laser  action  was  observed  by  Maiman  in  1960.     Many 
different  laser  systems  have  since  been  investigated.    One  of  the  most 
promising  systems  of  efficient  laser  energy  is  the  chemical  laser.     Since 
the  energy  is  stored  within  the  laser  medium  itself,  the  chemical  laser 
is  of  technological  interest  for  its  simplicity  and  compactness.     The 
ability  of  such  a  laser  to  provide  insight  into  the  distribution  of  chemical 
energy  in  the  products  of  a  reaction  make  it  an  attractive  tool  for  under- 
standing complex  reactions. 

The  most  favorable  reactions  for  the  production  of  laser  action 
are  those  which  are  highly  exothermic  and  involve  the  formation  of 
products  in  high  energy  states.     Such  reactions  are  most  commonly  of 
the  branched-chain  type.    The  excited  products  may  be  electronically 
or  vibrationally  excited.    However,  they  must  be  in  a  state  sufficiently 
differing  from  equilibrium  so  as  to  be  capable  of  producing  a  population 
inversion  either  within  the  products  or  by  energy  transfer  to  an  additive. 

The  latter  method  has  been  extensively  investigated  using  carbon 
dioxide  as  an  additive;  and  is  the  process  of  interest  in  the  Hydrogen 
Azide-Carbon  Dioxide  (HN     +  CO  )  laser.    The  explosive  decomposition 
of  HN     produces  vibrationally  excited  nitrogen  which  subsequently, 
through  near  resonant  energy  transfer,  produces  a  population  inversion 
within  the  vibrational  levels  of  CO    .    The  laser  action  then  occurs 
through  stimulated  emissions  from  the  asymmetric  vibration  level 
("V)  3  level)  of  C02  . 


In  the  following  sections  a  detailed  description  of  the  reaction 
mechanisms  involved  in  the  HN     +  CO     laser  will  be  presented.    The 

■J  Ld 

explosive  decomposition  of  HN     and  the  resulting  laser  gain  have  been 
analytically  modeled  to  gain  an  insight  into  the  potential  of  the  laser. 
In  addition,  an  experimental  investigation  has  been  conducted  in  order 
to  determine  the  existence  of  laser  action  and  applicability  of  the  model, 


II.     THEORY  OF  HYDROGEN  AZIDE  EXPLOSION 

Hydrogen  Azide  (HN   )  is  an  unstable  endothermic  gas  at  room 
temperature.     It  is  known  to  explode  violently  even  in  the  absence  of 
oxygen,  above  certain  limits  of  temperature  and  pressure.    At  low 
pressures,  less  than  40  mm  Hg ,  it  is  a  relatively  stable  compound. 
However,  the  stoichiometric  process  by  which  HN     reacts  at  these 
pressures  is  not  constant  (Ref.   1).    This  has  hampered  the  investigation 
and  understanding  of  the  mechanisms  by  which  the  reaction  proceeds. 

The  generally  adopted  reaction  mechanism  of  HN     is  (Refs .   1  -  7): 


HN3- 

— >    NH  +  N 

(a) 

N2 

+ 

HN3   - 

— ^  NH  +  2N 

(b) 

NH 

+ 

HN     - 
3 

— >  H      +  2N 
2             2 

(c) 

NH 

+ 

UN3- 

v.     -N1TT       _L     -Nl 

(d) 

>  NII2+  N3 

NH 

+ 

HN3- 

i.      -MfT      +    TU 

(dl) 

im2 

"  rJii3        3 

2N3    - 

->3N2 

(d2) 

NH 

+ 

NH   - 

— >H0       +  N0 

(e) 

Domination  of  reaction  (c)  results  in  a  violent  explosion  while  reaction 
(d)  leads  to  a  slowing  and  eventual  quenching  of  the  reaction.    The 
rates  at  which  the  two  mechanisms  proceed  are  a  function  of  the  pressure 
of  HN    ,  the  method  of  initiation,  and  the  type  and  quantity  of  any 
diluent  (Ref.   7)  . 
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Initiation  of  HN_  by  spark  discharge  and  flash  photolysis  lead  to 
a  domination  of  reaction  (c)  while  slow  thermal  decomposition  follows 
reaction  (d) .    However,  spark  discharge  and  flash  photolysis  initiation 
may  still  result  in  the  quenching  of  the  explosion  by  reaction  (d)  if  the 
dilution  of  HN     is  sufficiently  high. 

Several  investigators  have  suggested  reaction  (c)  be  divided  into 
reactions: 

NH  +HN3  — >    H2  +  2N2  (c) 

NH  +  HN     — »    H2N2  +  N2  — *  2NH  +  N£  (c') 

The  presence  of  H   N     has  not  been  experimentally  observed  and  reaction 
c'  is  5  Kcal  mole       less  exothermic  than  reaction  (c) .    In  addition, 
H   N_  is  a  relatively  stable  compound.     Therefore  it  is  assumed  that  if 
the  reaction  of  NH  with  HN     proceeds  via  a  highly  excited  H_N     particle, 
it  would  decompose  preferentially  into  H     +  N    . 

Reaction  (b)  as  a  propagation  step  has  not  been  observed  experi- 
mentally, however,  it  has  been  suggested  as  the  probable  mechanism 
(Ref  s .  2,3).    The  vibrational  excitation  of  nitrogen  is  consistent  with 
the  observed  lasing  action  obtained  from  the  reaction  of  HN     when 
diluted  with  CO     (Ref s  .   11,12).     It  is  energetically  feasible  since 
there  is  sufficient  excess  energy  from  flash  photolysis  and  reaction  (c) 
to  produce  vibrationally  excited  N     capable  of  reacting  with  HN 
(Ref.  2). 

V.  G.  Voronkov  et.  al.   (Ref.  8)  has  suggested  a  totally  different 
reaction  mechanism  in  which  the  propagation  and  branching  reactions 
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are  interchanged: 

NH  +  HN     — *-H     +  N   * 

N2*  +  fi  HN3  — *N2  +   y9HN  +  ^N2 

N    *  is  vibrationally  excited  and  capable  of  decomposing    fl  HN 
colecules.    These  reactions  are  plausible;  however,  A.  S.  Rozenberg 
and  V.  G.  Voronkov  (Ref s  .  9,  10)  in  later  publications  have  followed 
the  generally  accepted  reaction  mechanism.    Therefore  the  possibility 
of  the  interchange  of  these  reactions  has  not  been  considered. 
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III.    ANALYSIS  OF  HN     EXPLOSION 


If  the  concentration  limits  and  method  of  initiation  are  restricted, 
an  analysis  of  the  explosive  decomposition  of  HN     is  possible.     Since 
an  explosion  is  necessary  for  lasing  to  occur  and  the  method  of  initiation 
has  been  chosen  as  flash  photolysis,  prediction  of  the  products  of  the 
explosion  will  allow  the  laser  to  be  modeled. 

Provided  the  dilution  of  HN     is  not  sufficient  to  inhibit  an  explosion, 
reaction  (d)  may  be  neglected  since  the  products  of  the  reaction  have  not 
been  experimentally  observed  (Refs.    1,4).     The  equations  necessary  to 
describe  the  explosion  become: 

.HN3 >NH    +    N2*  (1) 

N    *    +  HN   >NH    +    N       +    N    '  (2) 

NH      +HN3 >H2     +    2N2*  (3) 

NH      +  NH    vH       +    N    '  (4) 

N    '     +  M    »N       +    M  (5) 

*%t    =    Htime.AE     )  (6) 

*  N     in  its  ninth  vibrational  level  or  higher 

1  N     in  less  than  its  ninth  vibrational  level 

NH  in  its  a' A      electronic  state 
These  reactions  are  divided  into  sections  according  to  their  number  so 
the  assumptions  and  calculations  necessary  to  the  model  may  be  stated 
in  detail. 
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(1)        HN     —  hS>  — >  NH  +  N    * 

Photochemical  decomposition  of  HN     was  first  investigated  by 

o 

Beckman  and  Dickerson  who  showed  the  quantum  yield  to  be  3.0  +  0.5 
and  independent  of  pressure  between  2  and  130  mm.  of  mercury.    Flash 
photolysis  of  HN     has  been  studied  by  Cornell,  Berry,  and  Lwowski 
(Ref .   13)  and  Thrush  (Ref .  4)  .    They  detected  the  presence  of  the  NH 
radical  and  identified  it  in  its  ground  electronic  state.     Okabe's  (Ref.  3) 
and  Welge's  (Ref s  .   14,   15)  recent  investigations  of  vacuum-uv  photolysis 

have  shown  the  primary  photochemical  process  in  the  dissociation  of  HN 

-I-  3      - 

to  result  in  N     and  NH  in  their  X'  3l      and  (a1  A     orX£"  )  states 
2  *— g 

respectively.    Okabe  also  reported  that  Stuhl  (Ref.   16)  detected  these 
as  the  primary  states  resulting  from  flash  photolysis.     Konar  et.  al. 
(Ref.   17)  has  suggested  N     to  be  in  a  vibrationally  excited  ground 
electronic  state  since  selection  rules  preclude  both  N     and  NH  from  being 
electronically  excited. 

It  is  assumed  that  NH  is  formed  in  its  a'A     state,  having  a  total 
energy  of  61  Kcal.  mole      .     For  an  average  quantum  wave  length  of 
2125°A  the  total  energy  available  for  dissociation  is  154  Kcal.  mole      , 
thus  leaving  an  average  of  93  Kcal.  mole       to  be  shared  between  the 
vibrational  states  of  NH  and  N?  .    The  principle  of  equal  parittion  of 
energy  indicates  that  nitrogen  should  be  in  its  ninth  vibrational  level, 
thus  having  an  excess  energy  of  60  Kcal.  mole      ,  more  than  sufficient 
for  decomposition  of  HN    .    The  remaining  13  Kcal.  are  assumed  to  be 
absorbed  as  translational  energy. 
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The  rate  constant  for  the  photolysis  of  HN     may  be  calculated  if 

the  absorption  coefficient,  quantum  yield  and  energy  of  the  flash  is 

known.     The  rate  of  decomposition  of  HN     depends  on  the  photon  flux 

o 

from  the  flashlamp,  its  absorption  coefficient  and  the  quantum  yield 
and  is  given  by  (Ref.    18): 


Rtjk      =    mv  tHN  J    =  0(1  -  S        )  Q 
HN3  3 

where  0  is  photon  flux,  Q  the  quantum  yield,  oc  the  absorption  coef- 


ficient, Jl   the  optical  thickness  of  the  absorber  and  KHV  the  average 
rate  constant. 

The  average  photon  flux  may  be  calculated  from: 

At  h\H3  A% 

where:  E     is  the  energy  stored  in  the  capacitor  bank 
o 

£  the  efficiency  of  the  flashlamp 

>^t  pulse  duration 

^   -fX  flash  spectral  range 


aO 
$ 


HN  absorption  band  of  HN 

HN  center  frequency  of  HN     absorption  band 

h  Planks  constant 

G  geometrical  relation  between  the  laser  tube  and 

flashlamp 

For  a  discharge  energy  of  lOOOj  the  ILC  7L24  lamp  has  an  efficiency  of 

85%  over  the  range  ^"^££  =  2000°A  to  4  microns.     The  absorption 

coefficient  of  HN     is  25  at  one  atmosphere  and  a  center  band  frequency 
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of  ""^HN     =  212  5°A.    The  absorption  band  (  A  "V  HN   )  of  HN     is 
between  2000°A  and  2250°A.    The  quantum  yield  is  2.0  (Ref.   5).    The 
geometrical  relation  between  the  flashlamp  and  laser  tube  is  based  on 
the  product  of  the  angle  which  the  laser  tube  subtends  on  the  flashlamp 
and  the  reflector  coupling  of  the  cavity  mirrors  .    The  angle  subtended 
by  the  laser  tube  utilized  in  this  inveatigation  was  0.2  of  the  total 

flash  angle  and  the  reflector  coupling  factor  was  assumed  to  be  2  . 5 

23 
(Ref.  19).    These  values  yield  a  photon  flux  of:  0  =  5.4  x  10       photons/sec 

The  average  rate  constant  for  the  decomposition  of  HN     may  be 

calculated  from  the  rate  equation: 


KHV  =  J2f  (1  -  e  )  Q/[HN  ] 

oCis  independent  of  pressure,  and  the  average  concentration  of  HN 

o 

1 8 
in  the  laser  tube  has  been  taken  as  10       particles/cc.     Substitution  of 

the  above  quantities  into  the  equation  yield  an  average  rate  constant  of 

2  -1 

3.23x10     sec 

From  the.  experimental  flash  profile  the  rate  of  decomposition  of 

HN     as  a  function  of  time  may  be  determined.    This  profile  is  shown  in 

Figure  15  and  may  best  be  described  by  a  cubic  equation: 

KHV  =  at3  +  bt2  +  ct  +  D 

From  the  average  rate  constant  and  the  experimental  profile  the  boundary 

conditions  on  the  above  equation  become: 
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t  =  0 

KHV  =  0 

t  =  100  usee .  , 

KHV  =  0 

t  —  37.5  usee  . 

r 

dKHV 
dt 

^X  =  100 

-)        KHVdt  = 

KHV 

tt  =  o 

The  above  boundary  conditions  yield  a  rate  constant  of 

15  3  11  2  7 

KHV  =2.89x10      t     -6.4x10     t     +3.42x10      t  where  t  is  in  seconds 

This  rate  constant  results  in  total  decomposition  of  HN     at  a  pressure  of 

6  torr.    When  2  5  torr  of  HN     that  is  sufficiently  diluted  to  prevent  an 

explosion,  only  9.5%  of  the  HN     is  decomposed. 

(2)        N*  +  HN3  — >NH  +  N2  +  N   ' 

The  reaction  of  N   *  and  HN     as  the  propagation  equation  for  the 
branched -chain  explosion  explains  the  observed  dependence  on  various 
diluents  and  hydrocarbon  scavengers  (Ref .  2) .     Inorganic  diluents  do 
not  significantly  reduce  the  rate  of  formation  of  products  for  volume 
ratios  of  diluent  to  HN     less  than  10.    For  ratios  greater  than  10  they 
slowly  decrease  the  rate  of  production  of  products  to  a  maximum  of  30% 
at  dilutions  of  100  and  greater.    Hydrocarbon  scavengers  however  have 
a  marked  effect  and  reduce  the  rate  by  70%  for  dilution  ratios  less  than 
10  (Ref.  2).     This  suggests  that  inorganic  diluents  do  not  react  with  the 
excited  NH  as  do  the  hydrocarbons.     Instead  they  merely  increase  the 
rate  of  relaxation  of  the  vibrationally  excited  N    .     Since  the  vibrational 
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levels  of  N„  below  the  ninth  are  not  sufficiently  excited  to  react  with 
HN    ,  this  increase  in  relaxation  rate  is  reflected  in  the  slowing  of 
reaction  (2) . 

Nitrogen  particles  in  vibrational  levels  above  the  ninth  have  been 
assumed  to  be  in  the  ninth  level.     The  number  of  particles  in  these 
levels  are  negligible  at  room  temperature.    At  higher  temperatures  the 
rate  of  relaxation  to  a  Boltzman  distribution  in  the  lower  levels  is  a 
minimum  of  ten  times  faster  than  the  excitation  to  levels  above  the 
ninth . 

The  rate  constant  for  the  reaction  of  HN     with  excited  nitrogen 

11        3  -1 

has  been  given  by  Okabe  (Ref.   3)  as  approximately  10       cm     moles 

sec       at  300°K.     In  order  to  obtain  the  temperature  dependence  of 

reaction  (2),  the  activation  energy  of  the  reaction  has  been  calculated 

— F  c\  /RT 
and  the  rate  is  assumed  to  obey  an  Arrhenius  equation  (K  =  Ae  ) . 

The  constant  A  may  then  be  determined  from  the  known  value  of  the  rate 

constant  at  300K. 

The  activation  energy  has  been  calculated  from  Simple  Collision 

Theory  (Ref.  20)  .    The  relative  Kinetic  energy  along  the  line  of  centers 

of  the  two  approaching  molecules  is  given  by: 

E 


1         r\,  2     2  2      2"] 

1    =    —  M     v       +  M.v. 

r        2ji      L      a      a  bb_J 


where:    M.  is  the  mass  of  reactant  i 
1 

V.  is  the  mean  molecular  velocity  of  i 
1 


JL 


is  the  reduced  mass 
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The  activation  energy  is  related  to  the  E    in  Simple  Collision  Theory  by: 

E     =  E    +    — 
a        r        2 

where:    T  =  temperature  in  °K 

R  =  universal  gas  constant 

The  mass  of  N     is  28.01  and  HN     is  43  .0.    At  a  temperature  of  300°K 

the  above  equations  yield  an  activation  energy  of  6.42  Kcal/mole.    The 

constant  A  can  then  be  obtained  from  the  Arrhenius  equation  at  a 

temperature  of  300°K: 

.       „  Ea/RT       7   c.  -9        3  -1  -1 

A  =  Ke  =  7.64  x  10       cm     part       sec 

and  the  rate  constant  can  then  be  expressed  as  a   function  of  temperature: 

„n       „    rr  -9     -6400/Rt        3  -1  -1 

K2  =  7.65  x  10       e  cm     part       sec 

(3)        NH  +  HN    • *H     +  2N* 

The  interaction  of  NH  with  HN     is  highly  exothermic  and  involves 
the  formation  of  several  bonds  which  is  indicative  of  a  branched-chain 
reaction,  and  the  creation  of  vibrationally  excited  particles  (Ref.  21). 
Since  the  energy  released  in  such  a  reaction  is  concentrated  along  the 
bonds  of  the  products,  and  with  the  previous  assumption  of  the  principle 
of  equal  partition  of  energy,  it  is  assumed  that  all  the  nitrogen  formed 
is  in  the  ninth  vibrational  level.    Hydrogen's  vibrational  levels  are  much 
more  widely  spaced  than  those  of  N     and  it  has  no  low-lying  electronic 
state  which  could  be  excited.    The  participation  of  hydrogen  as  a  re- 
actant  has  not  been  experimentally  observed  and  therefore  it  is  assumed 
that  if  hydrogen  is  formed  in  an  excited  state  its  energy  is  quickly 
dissipated  into  translational  energy. 
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The  rate  constant  for  the  reaction  of  HN     with  various  NH  radicals 

has  been  determined  by  Okabe  (Ref.3).    Assuming  all  NH  formed  to  be  in 

the  a' A     state  the  rate  of  reaction  (3)  is  5.5  x  10       cm     moles       sec 

at  300°K.    The  temperature  dependence  of  the  rate  constant  for  reaction 

(3)  has  been  determined  by  the  same  method  as  that  of  reaction  (2)  and 

has  the  value  of: 

K3  =  8.7x  10"7e     ^^  cc 

Ki 


part.  sec. 


(4)        NH+  NH =>H2  +  N2' 

Reaction  (4)  is  relatively  slow  and  may  be  neglected  in  the  pre- 
diction of  the  explosion  characteristic  of  HN      (Ref .    10)  .     It  does, 
however,  serve  to  explain  the  long  lasing  time  observed  from  an  ex- 
plosion of  HN     and  CO   ;  therefore  it  is  included  in  the  model.     The 
rate  constant  has  been  obtained  from  Reference  33  and  is: 

ita  o      c  in"12    m    °'55  -1900 

K4  =  3.6x10  T  e    -jj- 


part.  sec. 

(5)        N  #  +  M  »-N     +  M.(#=  any  vibrationally  excited  level) 

In  order  to  predict  the  effect  of  inorganic  diluents  on  the  explosion,  the 
collisional  deactivation  of  N  #  must  be  included.     Reaction  (5)  is  not 
kinetically  correct  but  merely  representative  of  the  process.     The  actual 
process  may  be  divided  in  two  parts:    collisional  relaxation  of  pure 
nitrogen  and  the  effects  of  additives. 

The  derivation  of  the  multi-quantum  level  population  equation 
is  shown  in  detail  in  Reference  22  and  will    not  be  repeated  here.    The 
only  assumption  involved  in  the  derivation  is  that  transitions  are  only 
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allowed  between  adjacent  quantum  levels.     The  validity  of  this  assumption 

is  supported  by  the  agreement  obtained  with  experimental  results  (Ref.  22). 

The  general  equation  for  the  vibrational  relaxation  of  pure  nitrogen  is 

given  by: 

dN 

-K  ,N     -K  nN    +K       ,       N         +K       ,       N 


dt  m,m+l    m     m,m-l    m     m+l,m    m+1      m-l,m    m-1 

where  K     is  the  transition  rate  for  N_  to  or  from  the  m      vibrational 
lj  2 

level. 

The  transition  rates  for  pure  nitrogen  have  been  obtained  from 
Reference  24  (translational-vibrational)  and  Reference  2  8  (vibrational- 
vibrational)  and  are: 

21    -1    (-290    4T~   '     -12    5) 
T-V  transitions  =  7.34x10     T      e^   ^u*qi  iZ,::);T=1950 

=7.34xl021T-1e(-56-7T"1/3-31'2)     T  =  1950 

1    1/2 

V-V  transitions  =  9.65x10   T 

where:    T  is  in  degrees  Kelvin  and  the  rates  in  sec 

The  increase  in  vibrational  relaxation  of  N     by  diluents  depends 
primarily  on  how  well  matched  the  vibrational  levels  of  the  diluents  are 
to  those  of  nitrogen.     The  process  of  vibrational  energy  transfer  for 
molecules  having  exactly  matched  levels  is  extremely  efficient,  occurring 
approximately  every  400  collisions,  while  the  probability  of  energy 
transfer  for  levels  differing  greater  than  200  cm       may  be  assumed 
negligible  for  most  reactions  (Ref.  22). 
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The    \?  7  level  of  carbon  dioxide  (CO   )  and  the     ^     level  of 
nitrogen  differ  by  only  18  cm       and  therefore  an  extremely  efficient 
near-resonant  energy  transfer  occurs  between  these  levels  .    The  effect 
on  the  explosion  by  CO     and  additional  diluents  of  interest  to  the  HN 
laser  are  shown  in  Figure  1. 

(6)        TEMPERATURE  AS  A  FUNCTION  OF  TIME 

The  final  quantity  necessary  for  prediction  of  the  products  of  the 
explosive  decomposition  of  HN     is  the  variation  of  the  translational 
temperature  with  time.    In  a  constant  volume  system  this  increase  in 
temperature  is  due  primarily  to  the  energy  evolved  from  the  reaction  into 
the  translational  and  rotational  motion.     Since  the  relaxation  of  the 
rotational  levels  of  a  molecule  are  characteristically  on  the  order  of 
nano-seconds,  the  energy  liberated  into  rotational  energy  may  be  con- 
sidered as  a  direct  increase  in  the  translational  temperature  (Ref s .  23,  24) 

The  equation  for  the  raise  in  translational  temperature  as  a  function 

of  time  may  then  be  expressed  as: 

5*  F  R 

d_T  ~    thermal   i 

dt  Z   Cv  [Xi] 

The  above  reactions  have  been  programmed  and  numerically  inte- 
grated on  the  IBM  360  computer.    The  resultant  temperature  from  the 
explosion  has  been  compared  with  the  experimental  results  of  A.  S. 
Rosenberg  and  V.  C.  Voronkov  (Ref.  9)  in  Figure  1.    Exact  agreement  is 
not  to  be  expected  since  the  experimental  data  were  obtained  from 
explosions  resulting  from  initiation  by  a  minimum  electrical  discharge. 


22 


IV.    ANALYSIS  OF  THE  HN     +  C02  LASER 


The  extension  of  the  HN     explosion  model  to  allow  calculation 
of  laser  gain  may  be  accomplished  relatively  easily.     The  primary 
equations  (effects  of  diluents  on  N   )  for  the  prediction  of  gain  are  con- 
tained in  the  explosion  model.     It  is  only  necessary  to  consider  in 
greater  detail  the  vibrational  relaxation  of  CO     and  its  energy  transfer 
with  the  vibrationally  excited  nitrogen  evolved  from  the  explosion. 

This  vibrational  energy  transfer  process  may  be  described  by 
dividing  the  vibrationally  excited  levels  of  CO     into  three  modes 
(Ref.  2  5).    As  shown  in  Figure  2,  mode  0  contains  only  CO     molecules 
in  their  ground  vibrational  level,  mode  I  those  in  the  symmetric  and 
bending  levels,  and  mode  II  being  those  in  the  asymmetric  stretching 
level. 

The  energy  of  the  asymmetric  stretching  level  (mode  II)  of  CO    , 
as  previously  noted,  is  near  resonant  with  the  v     level  of  N_  ,  therefore 
resulting  in  a  quantum  energy  transfer  approximately  every  465  collisions 
(Ref.  25).     In  addition,  N     in  its  ground  electronic  state  has  a  zero 
permanent  dipole  moment  and  cannot  decay  to  its  ground  vibrational 
level  through  electric-dipole  radiation.     Its  first  vibrational  level  is 
metastable,  thus  being  ideal  for  the  selective  excitation  of  the  asym- 
metric stretching  level  of  CO    . 

The  V     level  of  CO?  is  also  metastable,  having  (at  420°K  and  a 
total  pressure  of  15  torr)  a  radiative  lifetime  on  the  order  to  four  seconds 
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and  a  collisional  de-excition  rate  of  4  x  10       seconds  (Ref.  25).    The 

levels  of  mode  I  at  the  same  conditions  have  an  overall  relaxation  time 
of  approximately  2  0  ^-seconds,  being  limited  primarily  by  the  (Ol'O) 
level  (the  lowest  vibrational  level  of  CO_)  which  required  18    JU  -seconds 
for  deactivation.    This  large  difference  in  relaxation  times  allows  the 
inversion  mechanism  by  which  lasing  action  may  be  obtained. 

The  mechanism  of  stimulated  emission  (lasing  action)  in  the  CO 
molecule  is  characterized  by  transitions  between  two  vibrational  levels. 
The  upper  level  in  CO     is  the  asymmetric  stretching  level  (001)  while 
the  lower  level  is  the  symmetric  level  (100)  .     The  symmetry  of  carbon 
dioxide  allows  only  odd  rotational  levels  in  the  upper  vibrational  level 
and  only  even  in  the  lower  level  (Ref.  26).    Laser  transitions  occur 
between  particular  rotational  levels  within  the  vibrational  levels  ,  and 
correspond  to  a  rotational  quantum  change  of  J=+l .    Changes  of  J=+l 
are  P-branch  transitions  and  result  in  a  higher  gain  than  those  of  J=-l 
which  are  R-branch  transitions  .     Since  lasing  action  occurs  first  from 
the  rotational  level  of  the  highest  gain  which  is  subsequently  fed  by 

other  levels  to  maintain  a  Boltzman  distribution  (requiring  less  than 

-7 
10       sec),  R-branch  transitions  do  not  occur  (Ref.  27).    This  results 

in  a  laser  emission  at  a  wavelength  of  10.6  microns. 

In  order  to  calculate  the  energy  output  from  the  HN     +  CO     laser, 

the  number  of  molecules  in  each  vibrational  level  of  CO     must  be  known 

as  a  function  of  time.    However,  since  the  vibrational  levels  of  CO 

greater  than  (00°  1)  have  an  extremely  fast  relaxation  time  (200  n-sec  at 
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420°K  and  15  torr)  and  the  relaxation  of  the  lower  level  (100)  is  limited 

by  the  depopulation  of  the  (010)  level  to  the  ground  state,  the  population 

inversion  and  gain  may  be  calculated  within  5%  from  the  knowledge  of 

the  total  number  of  molecules  within  modes  I  and  II  (Ref.  2  5).     Considering 

the  assumptions  made  in  prediction  of  the  products  of  the  explosion  this 

is  sufficient  accuracy.    The  equations  for  the  population  inversion  and 

gain  are: 

PV  =  (N    A  ,     TT  -  N        ,     T)/N_^ 
mode  II         mode  I        C02 

Go         8TTKV2TTK      T3/2      A21  (2J+1)B(NII   NI)e  KT 

where:    PV  =  population  inversion 

G     =  gain 
o      3 

Other  symbols  have  their  usual  quantum  mechanical  meaning 
and  are  listed  in  the  TABLE  OF  SYMBOLS. 

The  reactions  and  rate  constants  which  must  be  added  to  the 
explosion  model  to  determine  the  above  quantities  are  given  in  Appendix 
A.    The  kinetic  rate  equations  have  been  numerically  integrated  using 
Fortran  IV  language  on  the  USNPGS  IBM  360  computer.    A  parameter 
variation  based  on  the  above  assumptions  has  been  carried  out  and  the 
results  plotted  in  Figures  3  and  4.    These  will  be  discussed  in  the 
Results  section. 
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V.     HN     +  CO     LASER  PROGRAM 


Equations  1-6  (Appendix  A)  result  in  a  set  of  first  order  linear 
differential  equations.     These  equations  may  be  integrated  to  determine 
the  products  as  a  function  of  time  if  the  rate  constants  and  initial  con- 
centrations of  the  reactants  are  specified.     This  integration  has  been 
carried  out  numerically  on  the  USNPGS  IBM  360  computer. 

The  program  and  a  sample  output  is  shown  in  Appendix  B.    The 
relation  between  computer  and  usual  notation  is  contained  on  comment 
cards  within  the  program.    A  region  of  100K  and  an  average  computer 
time  of  five  minutes  is  required  to  predict  laser  gain. 

The  USNPGS  RKLDEQ  library  subroutine,  employing  a  Runge- 
Kutta-Gill  predictor-corrector  method,  has  been  used  to  numerically 
integrate  the  rate  equations.    In  addition,  a  variable  step  size  based 
on  temperature  has  been  provided  to  reduce  computer  time.     Integration 
is  accomplished  on  a  particle  per  cubic  centimeter  basis  and  double 
precision  accuracy  (16  significant  digits)  has  been  employed  to  minimize 
mathematical  error. 

Data  output  is  controlled  by  numerically  recording  every  tenth 
step  while  laser  gain  and  population  inversion  are  graphically  displayed 
every  step  once  the  gain  becomes  greater  than  10      .    The  subroutine 
OSPLOT  has  been  used  for  plotting  gain  and  population  inversion.     It 
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allows  an  unlimited  number  of  points  to  be  plotted  and  requires  only  0.2 
seconds  of  computer  time  to  plot  an  average  of  2  00  points  . 

The  sample  program  in  Appendix  B  is  capable  of  calculating  laser 
gain,  population  inversion,  and  vibrational  temperature  for  the  HN     +  CO 

■J  £* 

laser  with  additional  diluents  of  helium,  nitrogen  or  hydrogen.     Lasing 
action  and  cavity  losses  have  not  been  considered  in  the  model.    If 
investigation  of  other  inorganic  diluents  is  desired,  it  may  be  accomplished 
by  setting  He  (Helium)  in  the  program  equal  to  the  diluents  concentration 
and  changing  the  rate  constants  specified  in  the  program,  to  those  of  the 
diluent  of  interest. 

The  results  obtained  (Appendix  B)  allow  determination  of  affects 
of  pressure  and  diluent  variation  on  the  laser  gain  and  population  inversion. 
The  assumptions  necessary  to  model  the  explosion  are  the  limiting  factors 
to  the  accuracy  of  the  results.    Experimental  determination  of  the  rate 
constants  for  reactions  2  and  3  would  resolve  this  limitation.    Experi- 
mental observation  of  laser  action  would  also  aid  in  determination  of  the 
quantitative  accuracy  of  the  model. 
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VI.     EXPERIMENTAL  TECHNIQUES 

A  chemical  laser  requires  the  integration  of  several  very  different 
systems.    These  apparatus  may  be  divided  into  five  categories: 
(a)  general  description;   (b)  flash  equipment;   (c)  chemical  and  vacuum 
systems;   (d)  laser  cavity;  (e)  recording  equipment.    This  allows  a  more 
detailed  description  of  each  system  and  its  function  in  a  chemical  laser. 

A.         GENERAL  DESCRIPTION 

A  general  block  diagram  of  the  laser  components  is  shown  in 
Figure  5.    The  flashlamp,  laser  cavity  and  toxic  gas  systems  are  located 
in  a  fume  hood  due  to  the  toxic  and  explosive  nature  of  the  chemicals 
involved  in  the  generation  of  HN    .    The  fume  hood  is  exhausted  to  the 
atmosphere  by  a  fan  and  has  been  fitted  with  a  plexiglass  covered 
extension  for  greater  ease  in  weighing  chemicals. 

To  the  left  of  the  fume  hood  is  the  main  control  panel  and  the 
remainder  of  the  gas  system  containing  non-toxic,  non-reactive  gases. 
All  systems  can  be  operated  from  the  main  panel  with  the  exception  of 
charging  the  circuit  capacitor  bank.     In  addition,  all  recording  equip- 
ment is  located  in  this  area  for  greater  safety.     On  the  right  is  the 
trigger  circuit,  electrical  control  box,  capacitor  bank,  and  the  charging 
circuit  for  the  flash  system. 

Considering  the  high  voltages  involved  in  the  flash  equipment 
and  the  explosive  nature  of  the  chemicals,  the  two  systems  have  been 
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separated  wherever  possible.     During  all  experiments  the  operator  is 
located  in  area  to  the  left  of  the  fume  hood  to  minimize  contact  with 
the  high  voltage  systems  .    The  extension  to  the  fume  hood  is  fitted  with 
rubber  seals  and  in  combination  with  the  exhaust  fan  provides  a  sub- 
atmospheric  pressure  in  the  interior  insuring  an  air  flow  away  from  the 
operator.    The  safety  features  of  each  system  will  be  covered  in  the 
detailed  description  of  the  categories  below. 

B.         FLASH  EQUIPMENT 

The  flash  system  contains  three  electrically  separated  circuits: 
the  discharge  circuit,  triggering  circuit,  and  charging  circuit.    The 
three  components  are  integrated  through  the  electrical  control  box 
(Fig.  6)  located  to  the  right  of  the  fume  hood.    The  entire  flash  system 
can  be  monitored  and  fired  through  the  control  box.     However,  an 
additional  firing  switch  has  been  installed  on  the  main  control  panel 
for  greater  safety. 

The  discharge  circuit  (Fig.   7)  is  a  simple  series  circuit  containing 
the  flashlamp,  capacitor  bank  and  inductance  coil.     These  may  be  varied 
in  accordance  with  Reference  30  to  obtain  the  desired  pulse  duration 
and  energy  output.    The  flashlamp  parameters  specified  by  the  manu- 
facturer determine  the  bounds  within  which  the  system  may  be  operated. 
The  flashlamp  is  a  xenon  filled,  ILC  model  7L24  (60  cm.  by  7  mm.), 
lamp  with  nickel  plated  copper  electrodes.    A  lamp  with  either  a  natural 
or  synthetic  quartz  envelope  may  be  used  in  the  circuit  as  designed. 
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The  former  emits  light  of  2000°A  and  greater  while  the  latter's  emission 
band  starts  at  1600°A.    The  synthetic  quartz  lamp  is  particularly  useful 
for  decomposition  of  HN     which  has  an  absorption  band  that  rises 
sharply  below  2000°A.    However,  in  order  to  obtain  the  full  effectiveness 
from  the  lamp  the  laser  tube  must  also  be  made  of  synthetic  quartz.    The 
7L24  lamp  with  either  envelope  has  an  impedence  factor  of  116  and 
requires  a  minimum  trigger  voltage  of  25KV  for  2.4  micro-seconds. 

The  other. two  components  of  the  system  are  the  capacitor  bank, 
consisting  of  four  7  /Li.  and  two  1.5  yUf.  high  voltage  (2  5KV  capability) 
capacitors  in  parallel,  and  a  30.0  Jbcb.  inductance  coil  made  of  3/8-inch 
copper  tubing  consisting  of  25  windings  15  cm.  in  diameter.    The 
components  have  been  silver  soldered  in  series  using  3  x  1/4-inch 
copper  bars  and  3/8-inch  copper  tubing.     The  measured  ohmic  resistance 
of  the  circuit  is  0.0013  ohms  and  can  be  neglected  in  the  calculation  of 
the  operating  parameters. 

The  flash  circuit  has  been  designed  for  a  damping  factor  of  0.8 
and  a  pulse  duration  of  100;u-sec,  with  an  energy  output  of  2  000  joules. 
The  lamp  may  be  operated  under  manufacturer's  warranty,  with  a  damping 
factor  of  1.1  to  0.7.     This  allows  a  variation  of  energy  output  with  a 
minimum  increase  or  decrease  in  pulse  duration.     The  energy  output 
variation  with  capacitor  bank  voltage  has  been  calculated  from  Refer- 
ence 30  for  the  allowable  damping  factor  range.    Figure  8  shows  the 
energy  output  in  joules  for  a  given  capacitor  bank  voltage. 
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The  triggering  circuit  creates  an  ionized  spark  streamer  between 
the  two  electrodes  which  allows  the  main  discharge  to  occur.     It  is 
external  to  the  flashlamp  (Fig.  9)  and  consists  of  a  Dressen-Branes 
300  volt,   70    jla  DC  power  supply,  an  8   jjlI.  oil  filled  capacitor, 
SCR,  and  a  ILC  model  T105  trigger  transformer  with  a  step-up  ratio  of 
60:1.    The  circuit  is  controlled  by  a  relay  that  is  energized  by  pressing 
the  firing  switch.    A  detailed  circuit  diagram  is  shown  in  Figure  9. 

The  charging  circuit  consists  of  a  NJE  (New  Jersey  Electronics) 
Model  HA-51  variable  high  voltage  power  supply  (0-30kv,   0-lQjUa)  which 
is  connected  through  a  high  voltage  vacuum  switch  to  the  capacitor  bank 
in  series  with  a  3mjX     resistor  (Fig.   10).    When  the  charge/dump  switch 
is  placed  in  the  charge  position,  the  capacitor  bank  voltage  is  monitored 
on  a  20    ^aftAPI  (Assemble  Products,  Inc.)  meter  relay  on  the  electrical 
control  box.    The  meter  is  calibrated  as  a  voltmeter  and  interrupts  the 
charging,  by  isolating  the  charging  circuit  and  placing  the  system  in 
the  ready-fire  position,  when  the  set  voltage  is  reached.    A  warning 
flasher  is  mounted  on  the  control  box  and  is  activated  at  any  time  the 
charge/dump  switch  is  placed  in  the  charge  position. 

Safety  is  of  prime  importance  to  the  flash  system.    With  the 
exception  of  the  flashlamp,  the  entire  system  has  been  placed  on  the 
side  of  the  laboratory  away  from  the  operator's  position  during  experi- 
ments.    The  flasher  on  the  electrical  control  box  is  automatically 
activated  when  charging  is  initiated.    The  relay  meter  which  separates 
the  charging  system  when  the  desired  voltage  is  reached  must  be 
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electrically  unslaved  after  firing  before  the  system  may  be  recharged. 
Plywood  covers  have  been  placed  over  the  capacitor  bank,  inductance 
coil  and  high-voltage  switches  to  minimize  the  danger  of  operator  contact 
or  possible  arcing.    As  constructed,  the  system  is  easily  operable  while 
providing  maximum  operator  safety. 

C.         CHEMICAL  AND  VACUUM  SYSTEMS 

The  chemical  and  vacuum  systems  have  been  integrated  through 
a  manifold  mounted  directly  above  the  laser  cavity  in  the  fume  hood 
(Fig.   11).    All  components  of  the  toxic  gas  generator  are  contained  under 
the  fume  hood  while  the  inert  gases  and  vacuum  pump  are  external  to 
the  fume  hood  and  connected  to  the  manifold  by  1/8-inch  stainless 
steel  tubing.    The  minifold  consists  of  stainless  steel  valves  that  are 
operated  by  solenoids  from  the  main  control  panel.    This  allows  a  fast 
and  easy  way  to  isolate  any  part  of  the  system. 

The  toxic/explosive  gas  generator  has  been  designed  primarily 
for  use  in  the  generation  of  Hydrozoic  Acid  (HN   )  .    However,  it  is 
capable  of  being  used  to  generate  other  chemicals  in  a  gaseous  state. 
The  description  here  will  be  concerned  with  the  generation  of  HN     since 
its  design  was  for  this  purpose. 

The  gas  generator  is  a  500  ml.  pyrex  flask  (Fig.   12)  fitted  with 
a  pyrex  berrette  graduated  in  cubic  millimeters  .    A  relief  valve  has  been 
provided  which  opens  to  exhaust  at  atmospheric  pressure  to  prevent 
possible  damage  from  over  pressure  during  generation  of  HN    .     The 
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system  is  connected  to  the  manifold  through  a  1/4 -inch  pyrex  U-tube 
by  a  combination  of  two  ball  joints  and  a  flexible  connection  (bellows) . 
The  U-tube  provides  a  convenient  way  of  drying  the  HN     with  phos- 
phorus  pentoxide  before  it  is  introduced  into  the  laser  tube. 

The  generation  of  HN     is  accomplished  by  the  following  reaction 

O 

(Ref.  31): 

NaN     +H   P04— »HN     f    + NaH   PO 

o  o  v3  &*  ^ 

Sodium  Azide  (NaN    :    powdered  solid)  and  phosphoric  acid  (H   PO   : 
85%  syrupy  solution  in  water)  are  mixed  according  to  their  gram  atomic 
weights  to  provide  the  desired  amount  or  pressure  of  HN     for  the  volume 

O 

of  the  generation  system.    Table  I  shows  the  weights  necessary  to 
generate  HN     to  pressures  of  50mm,  based  on  a  one-liter  system.     If 

O 

the  present  one-liter  system  is  changed,  the  quantities  of  NaN     and 

H   PO.  necessary  for  generating  the  desired  pressure  of  HN,,  may  be 

obtained  by  multiplying  the  data  in  Table  I  by  the  volume  of  the  system. 

To  generate  a  larger  quantity  of  HN     at  the  same  pressure,  the  laser 

tube  and  pressure  gauge  may  be  left  open  during  generation.    This 

increases  the  volume  from  one  to  six  and  one-half  liters. 

Sodium  azide  is  weighed  and  placed  in  the  generation  flask 

through  the  relief  valve  assembly  and  then  mixed  with  Dow  Corning 

704  diffusion  pump  fluid  (silicone  based)  DCC  (Dow  Corning  Corp.)  by 

means  of  a  teflon  coated  magnetic  stirrer  in  the  bottom  of  the  flask.    The 

quantity  of  H   PO.  used  in  the  generation  of  HN     is  chosen  slightly 

higher  than  stoichiometric  to  insure  all  NaN      (explosive  and  toxic) 

o 
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is  consumed.    After  the  system  has  been  evacuated  the  H   PO     is  intro- 
duced  into  the  generation  flask,  by  a  mechanically  operated  stopcock. 
The  stopcock  is  controlled  from  a  switch  on  the  main  control  panel  and 
rotates  from  fully  open  to  fully  closed  in  15  seconds. 

A  plexiglass  covered  extension  to  the  fume  hood  has  been  con- 
structed to  allow  greater  safety  in  the  weighing  and  handling  of  the 
chemicals  (Fig.  5) .    The  extension  has  been  sealed  by  the  use  of  rubber 
stripping  along  the  edges.    A  14-  by  12-inch  door  on  one  side  allows 
the  operator,  wearing  a  face  shield,  easy  access  to  all  parts  of  the 
system  while  handling  the  chemicals.     The  exhaust  fan  provides  an  air- 
flow away  from  the  operator,  which  is  exhausted  to  the  atmosphere,  to 

prevent  inhalation  of  toxic  HN     fumes. 

o 

The  inert  gas  system  is  external  to  the  fume  hood  and  is  connected 
to  the  manifold  by  1/8-inch  stainless  steel  tubing.    The  gases  used  in 
the  laser  were  manufactured  by  ARC  (Air  Reduction  Co.  ,  Inc.)  and  are 
supplied  in  one-liter  pyrex  flasks  at  a  pressure  of  755  mm  of  mercury. 
Their  guaranteed  purity  is: 

Carbon  Dioxide  CO        :99.995% 

Nitrogen  N  :99.999% 

Helium  He         :99.9995% 

Each  gas  is  introduced  into  the  laser  tube  by  a  separate  solenoid  switch. 
The  quantity  and  rate  of  induction  is  controlled  by  needle  valves  located 
immediately  behind  the  inlets.    The  pressure  in  the  laser  tube  is  moni- 
tored by  a  stainesss  steel  bourdon  tube  guage  with  a  range  from  0  to 

800  mm  of  mercury  and  a  precision  of  0.1%  full  scale. 
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Prior  to  introducing  the  laser  gas  mixture  the  system  is  evacuated 

_3 
by  a  mechanical  vacuum  pump  capable  of  producing  a  vacuum  of  10 

torr.    The  system  is  flushed  several  times  during  evacuation  with  an 

inert  gas  (N   )  to  insure  a  minimum  of  impurities  in  the  laser  tube.    Once 

evacuation  is  complete  the  desired  mixture  of  HN     and  additional  gases 

are  introduced  into  the  laser  tube. 

D.         LASER  AND  OPTICAL  CAVITY 

The  laser  components  are  mounted  on  a  one-meter  standard 
optical  bench  (Figs.   13,    14).    The  laser  tube  is  connected  to  the  mani- 
fold, via  a  ball  and  socket  joint,  by  1/4-inch  glass  tubing  which  is 
sealed  to  a  stainless  steel  bellows. 

The  laser  tube  is  made  of  fused  quartz  and  is  60  cm.  long  with  an 
inner  diameter  of  30. mm.    The  ends  of  the  tube  have  been  cut  and  ground 
at  a  Brewster  angle  of  56°  40'   (for  sodium  chloride)  to  minimize  reflec- 
tions.   Sodium  chloride  windows  were  attached  to  the  ends  of  the  quartz 
tube  by  chemically  inert  silicone  based  cement. 

The  optical  cavity  is  formed  by  two  gold-surfaced  concaved 
mirrors  one  inch  in  diameter  placed  78  cm.  apart  in  gimbal  mounts  made 
by  Oleo  Optical  Co.    The  mirrors  have  a  radius  of  curvature  of  ten 
meters.    A  2 -mm.  hole  drilled  in  the  center  of  one  mirror  provides  an 
output  for  the  recording  equipment  and  allows  easy  access  to  the  cavity 
for  alignment. 
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The  amount  of  energy  removed  by  the  2-mm.  hole  may  be  calcu- 
lated from  Gaussion  theory  and  is  directly  related  to  the  laser  beam 
spot  size  (Ref.   32).    This  becomes  important  when  small  gains  are 
involved  and  will  be  discussed  in  the  results  section.    For  the  laser 
cavity  described  above,  Gaussion  theory  yields  a  spot  size  6.0  mm.  in 
diameter  for  a  wavelength  of  10.6  microns  and  energy  removal  of  37%. 
This  applies  to  the  lowest  order  mode.     Higher  order  modes  which  exist 
in  the  comparatively  large  diameter  laser  tube  extend  further  in  the 
transverse  direction,  thus  most  of  their  energy  is  off  the  axis  of  the 
resonator. 

Alignment  of  the  laser  cavity  and  recording  equipment  is  compar- 
atively critical  and  conditions  for  stability  of  the  cavity  are  given  below: 

o  =  8,  g2  s  i 

L  =    RlR2 

Where 

i        L 
gl  =  l  ~    Rx 

i         L 

L  =  length  of  the  cavity 

R]  =  radius  of  the  mirror 
The  alignment  is  accomplished  using  a  commercial  He-Ne  laser  with 
a  continuous  output  of  one  milliwatt.    The  He-Ne  laser  is  placed  in  line 
with  the  laser  cavity  and  its  beam  is  passed  through  the  2  mm.  hole  in 
one  of  the  cavities'  mirrors.    The  mirrors  are  then  adjusted  until  all 
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reflections  within  the  cavity  pass  back  through  the  hole  and  directly 

-4 
into  the  He-Ne  laser.    The  limit  to  the  precision  of  alignment  (2  x  10 

radians)  is  the  smallest  adjustment  which  can  be  made  with  the  gimbal's 

micrometers . 

E.         RECORDING  EQUIPMENT 

In  order  to  effectively  monitor  the  light  output  from  the  laser 
cavity  two  detectors  were  used.    This  allowed  the  flashlamp  discharge 
to  be  recorded  by  one  detector  while  the  other  was  being  used  in  con- 
junction with  various  filters  to  determine  any  resulting  laser  energy. 
The  two  detector  signals  were  recorded  on  a  dual  beam  oscilloscope 
during  each  experiment.     The  position  of  the  two  detectors  in  the  experi- 
mental setup  is  shown  in  Figures  13  and  14. 

A  photo-diode  detector  was  used  to  monitor  the  flashlamp  discharge. 
It  provided  an  effective  method  of  monitoring  the  lamp  pulse  duration 
and  relative  magnitude  for  various  capacitor  bank  voltages.     In  addition, 
it  served  as  the  triggering  device  for  the  oscilloscope,  thus  providing 
a  reference  for  the  second  detector. 

For  detection  of  laser  energy  a  gold-doped  germanium  detector 
was  used  in  conjunction  with  several  filters.     The  detector  is  cooled 
by  liquid  nitrogen  and  biased  controlled  to  yield  a  signal  of  several 
millivolts.     Its  maximum  response  time  is  50  nanoseconds  over 
a  range  of  1.0  to  10.0  microns.     Because  of  the  detectors  high  sensitivity, 
a  monochromator  and  several  different  filters  were  employed  in  order  to 
eliminate  extraneous  light. 
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A    Perkrin-Elmer  Grating  monochromator  (MG-12b)  with  a  blaze 
angle  of  12  microns  was  used  to  monitor  a  specific  wavelength  of  light. 
It  was  calibrated  with  a  He-Ne  laser  and  cross-checked  with  a  mercury 
light.    When  the  doors  of  the  monochromator  are  fully  open  (2  mm.),  a 
wavelength  band  of  0  . 3  microns  was  provided  to  the  detector,  thus 
eliminating  the  requirement  of  filters  to  attenuate  the  light  signal. 

Two  different  filters  were  employed  without  the  monochromator  in 
an  effort  to  determine  lasing  action  that  might  have  occurred  in  a  range 
not  allowed  by  the  monochromator.    First,  an  Irtran  II  filter  which  passes 
light  between  2.0  and  12  microns  was  used.     However,  due  to  the  amount 
of  light  emitted  in  this  range  by  the  flashlamp,  the  oscilloscope  vertical 
scale  deflection  was  limited  to  50  millivolts  per  cm.     This  was  an  order 
of  magnitude  greater  than  the  expect  laser  signal  and  limited  determination 
of  any  output.    A  KBR  filter  passing  light  only  above  10.2  microns  was 
also  employed.     It  allowed  oscilloscope  deflection  as  small  as  5 
millivolts  per  cm.  but  eliminated  detection  of  the  9.4  micron  band  of 

co2. 

A  sample  output  of  each  recording  setup  is  shown  in  Figures  15 
and  16.    The  oscilloscope  vertical  response  time  is  0 .  06  ^-seconds 
and  scale  deflection  settings  are  shown  beside  each  trace.    The  minimum 
detectable  response  is  considered  to  be  five  millivolts  in  all  cases,  due 
to  radiation  noise  from  the  flash  circuit  discharge. 
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VII.    RESULTS 

A.         THEORETICAL 

The  agreement  of  the  computer  model  with  the  experimental 
adiabatic  temperature  (Fig.    1)  obtained  in  Reference  9  showed  maximum 
difference  of  9%.    This  lends  some  support  to  the  many  assumptions 
necessary  to  model  the  explosion.    Since  the  rate  constants  for  reactions 
(b  and  c)  had  to.be  calculated  from  sparse  experimental  data,  the  deter- 
mination of  exact  rate  constant  as  a  function  of  temperature  would 
considerably  strengthen  the  model. 

From  the  large  amount  of  data  obtained  from  the  computer  model, 
three  factors  are  of  major  interest  to  the  experimental  investigation  of 
the  laser.    First  is  the  magnitude  of  the  laser  gain  for  mixtures  of 
HN     +  CO-.     Second,  the  duration  for  which  the  gain  is  sufficiently 
high  that  it  may  be  experimentally  observed.    And  last,  is  the  effect 
of  the  addition  of  He  to  the  laser  mixture.    These  factors  have  been 
graphically  displayed  in  Figures  3  and  4. 

Figure  3  shows  laser  gain  as  a  function  of  the  dilution  ratio  of 
CO     to  HN    .    The  ratios  for  which  gain  is  obtainable  is  broadened  as 
the  pressure  of  HN„  is  reduced.    This  is  indicative  of  a  total  pressure 
limitation  on  the  laser.    For  40  torr  of  HN     and  a  total  pressure  of 
150  torr  the  vibrational  relaxation  rate  of  CO     becomes  comparable  to 
the  reaction  rates  of  the  explosion  which  precludes  laser  gain.     In 
addition,  the  maximum  pulse  duration  occurs  at  dilution  ratios  less 
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than  those  of  maximum  gain.    The  pulse  duration  falls  slowly  on  either 
side  of  the  maximum.    This  becomes  important  in  the  experimental 
observation  of  laser  action  only  at  high  pressures  where  the  pulse 
duration  approaches  the  response  time  of  the  experimental  recording 
equipment. 

The  effects  of  the  addition  of  He  to  the  laser  mixture  are  quite 
pronounced.    A  significant  increase  in  gain  as  well  as  pulse  duration 
results.    Figure  4  shows  increase  in  gain  with  dilution  of  He  for  25  torr 
of  HN    .    With  100  torrs  He  added  to  the  mixture,  the  gain  is  more  than 
doubled  and  the  maximum  pulse  increases  from  5  to  2  6^60.    Further 
increases  in  He  results  in  a  slight  decrease  in  gain  while  the  pulse 
duration  remains  nearly  constant.     However,  dilutions  of  greater  than 
4  to  1  for  He  is  approaching  the  lower  limit  of  the  model  where  the 
reactions  neglected  become  important. 

From  the  results  of  the  computer  model,  it  would  appear  that  laser 

action  could  be  obtained  even  in  the  crudest  of  experimental  setups. 

However,  two  things  have  not  been  considered  in  modeling  the  HN     + 

CO     laser.    The  actual  particle  depletion  that  results  from  laser  radiation 

and  the  possibility  of  a  reaction  involving  NH  and  CO    .    The  addition  of 

the  former  requires  a  significant  increase  in  the  number  of  assumptions 

necessary  for  modeling  the  laser  and  could  well  degrade  the  accuracy 

of  the  results  .    The  latter  would  only  be  speculation  since  very  little 

is  known  about  the  reaction  of  NH  with  any  molecule.    The  limitations 

and  accuracy  of  the  model  are  contingent  on  experimental  observation 

of  laser  action. 
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B.         EXPERIMENTAL 

An  experimental  investigation  of  the  HN~  +  CO?  laser  was  con- 
ducted for  pressures  of  HN     from  6  to  35  torr.    In  addition  to  pure  C09  , 
Helium  (He)  was  also  used  as  a  diluent  under  CO     starved  conditions 
in  order  to  maintain  a  temperature  less  than  2500°K.    Although  no  laser 
action  was  observed,  it  was  confirmed  that  flash  photolysis  in  excess 
of  1000J.  produced  explosive  decomposition  of  HN     and  resulted  in  a 
100%  increase  in  pressure  from  HN    .    This  result  adds  some  support  to 
the  assumptions  of  the  explosion  model  but  does  not  provide  an  insight 
into  the  reasons  for  the  unobserved  laser  action. 

Of  primary  importance  to  laser  action  is  the  alignment  of  the 
optical  cavity.     The  alignment  of  the  far  cavity  mirror  is  certainly  within 
the  limits  necessary;    however,  the  relatively  large  diameter  hole 
(2.0  mm.)  in  the  output  mirror  makes  its  alighment  difficult.     In  addition 
the  vibrations  from  the  vacuum  pump,  which  was  mounted  near  the  laser 
cavity,  adds  to  the  doubt  of  the  necessary  accuracy  in  alignment. 

The  large  amount  of  energy  (37%)  extracted  by  a  2 -mm.  diameter 
hole,  in  addition  to  the  Brewster  window  and  mirror  losses   (estimated 
at  2.0%),  would  preclude  laser  action  if  the  gain  were  lower  than  that 
predicted  by  the  model.    For  this  reason,  two  closed  mirrors  were  also 
used;  however,  alignment  was  even  more  doubtful  under  this  setup  since 
it  required  alignment  on  the  front  of  one  mirror  and  the  back  of  the  other. 


41 


Two  additional  factors  which  hampered  the  detection  of  laser 
action  were  radiation  from  the  triggering  circuit  and  the  necessity  of 
storing  the  laser  mixture  in  the  system's  only  pressure  gage.    The  radi- 
ation was  on  the  order  of  5  millivolts  per  centimeter  and  therefore  would 
only  be  significant  under  conditions  of  small  gain.    The  trigger  circuit 
was  redesigned  to  no  avail  and  finally  a  pulse  generator  was  ordered; 
however,  delays  in  delivery  prevented  its  use.    The  storing  of  the  laser 
mixture  in  the  pressure  gage,  however,  precluded  its  use  during  evacu- 
ation of  the  laser  tube  after  each  experiment.    The  elimination  of  the 
products  from  the  laser  tube  therefore  was  based  on  time.     However, 
this  is  thought  to  be  of  minor  importance  and  is  easily  remedied  by  the 
addition  of  a  small  pressure  sensing  device. 

The  latter  problems  are  of  minor  importance  and  correction  was  only 
limited  by  the  time  required  for  delivery  of  the  necessary  equipment. 
However,  alignment  is  most  critical  to  the  detection  of  laser  action. 
For  this  reason  it  is  thought  that  a  beam  splitter  internal  to  the  resonator 
cavity  should  be  used  for  alignment  so  both  mirrors  may  be  closed  and 
aligned  at  the  same  time.    The  most  effective  method  of  determination 
of  laser  output  is  that  of  a  gain  experiment,  and  in  addition  it  insures 
alignment  of  all  components  in  the  system  at  the  same  time.     The  timely 
repair  of  the  CO     laser  which  was  to  be  used  in  gain  experimentation 
prevented  this  type  of  investigation. 
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VIII.     CONCLUSIONS  AND  RECOMMENDATIONS 

The  major  objectives  of  the  research  conducted  were  the  analytical 
modeling  and  experimental  determination  of  laser  action  for  the  HN     +  C00 
chemical  laser.    The  analytical  model  has  been  completed  in  the  form  of 
a  computer  model  which  predicts  the  magnitude  and  duration  of  laser 
gain  as  well  as  the  population  inversion  within  the  laser  medium.    Experi- 
mental observation  of  laser  action  has  not  been  accomplished;  however, 
the  results  obtained  indicate  difficulties  in  the  experimental  setup  to  be 
the  reason. 

The  following  recommendations  are  made  for  possible  improvements 
in  the  model  and  experimental  investigation: 

A.  ANALYTICAL  MODEL 

1.  Experimental  determination  of  the  rate  constants  for 
reactions  b  and  c. 

2.  Inclusion  of  the  reactions  resulting  in  the  formation 
of  NH    . 

3.  Expansion  of  the  model  to  include  lasing  energy 
produced  and  percentage  of  energy  lost  due  to  the 
experimental  setup. 

B.  EXPERIMENTAL 

1.         A  gain  experiment  to  determine  the  existence  of 
laser  action. 
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2.  The  use  of  a  partial  transmitting  mirror  as  part  of  the 
laser  cavity  to  reduce  the  amount  of  energy  removed 
from  the  cavity  on  each  pass  . 

3.  Procurement  of  a  Trigger  Pulse  Generator  to  reduce 
the  electrical  noise  picked  up  by  the  recording 
equipment. 

4.  Inclusion  of  a  storage  system  for  HN     in  the  experi- 
mental  setup.    This  would  reduce  exposure  to  the 
toxic  chemicals  used  to  generate  HN     to  allow  more 
experiments  per  generation. 

5.  Extension  of  the  laser  cavity  to  increase  the  laser 
gain  on  each  pass  through  the  cavity. 

6.  The  use  of  a  rapid  frequency  scanner  to  determine  the 
wavelength  of  laser  action. 

7.  Operation  of  the  laser  tube  as  a  pure  CO     +  N     laser 
by  high  voltage  excitation  would  also  provide  the 
necessary  alignment  information. 
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U-tube 


Figure  12.  HN~  Gas  Generator 
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50  usec/cm 
Figurel5»  Oscilloscope  trace  of  output  from  Photo  Diode 
(upper)  and  IR  Detector  with  10.2  filter  (lower). 


50  usec/cm 
Figure  Ifc.  Oscilloscope  trace  of  output  from  Photo  Diode 
(upper)  and  IR  Detector  with  Irtran  II  filter  (lower). 


58 


HN3  torr 

NaN3  grams 

H3F0^  grams 

10 

0.0332 

0.0576 

20 

0.076^ 

0.1153 

30 

0.1132 

0.1729 

^0 

0.1520 

0.2305 

50 

0.1900 

0.28S0 

TABLE  I.   Amounts  of  NaN?  and 


H~P0^  necessary  to  generate  a 
desired  pressure 
one'  liter  system. 


desired  pressure  of  }ti-!~   in  a 
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APPENDIX  A:    REACTIONS  AND  RATE  CONSTANTS  USED  IN  THE  COMPUTER 
MODEL 


All  rate  constants  are  in  the  units  of:  part,    cc.sec 

HN     *  NH  +  N* 

3  2 

KHV  =  2.89  x  1015t3  -6.4  x  lC^V  +  3.42  x  10?t 

N*  +  HN     — *  NH  +  N^    +  N 

v0  .  -9    -6400/RT 

K2  =  7  .  64  x  10      e 

NH  +  HN     — »  H     +  2N* 

O  Cd  Li 

vi  -  o    7       in-7    "8210/RT 
K3  =  8 .  7  x  10      e 

NH  +  NH  — »  H     +  N' 

K4  =  3.6xl0-12T°-5V1900/RT  (Ref.  33) 

it 

N^  +  M  —>  N     +  M 

This  equation  as  stated  before  is  only  a  representation  of  the 
reactions  that  must  be  considered  in  the  computer  model.    The  following 
equations  are  those  actually  required  for  each  diluent  considered. 

N'2  +  H2  "^  N2  +  H2 

-1/3 
-8.76T     '    -2.23 
K  =  e 

N     (MODE  II)  +  CO     «--  *  N     +  CO    (MODE  II) 


,1/2 
K  =  e  -..--       forT^1000°K 


-0.0472T   '     -27.36 


1/2 
K  =  e0.l53T        -33.83  for  T  ^  1000°K  (Ref.  29) 
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N    (MODE  II)  +  He  --»  N     +  He 


K=1.36xl0-16T10-(60-7T"1/3-4-168)  (Ref.  24) 


C02(MODE  I)  +  C02  -->  C02  +  C02 

K^.SexlO-^TlO^17-81"173-1'808)  (Ref.  24) 


C02(MODE  I)  +  H2  — *  C02  +  H2 

K  =  4.0x10"  (Ref.  29) 


C02(MODE  I)  +  N2  — *  C02  +  N2 

K^^xlO-^TlO^17-81"173"1-80^  (Ref.  24) 


CO    (MODE  II)  +  He  — >  CO     +  He 

K  =  3-46.23T-1/3-23.01  (Rrf_ 


CO    (MODE  II)  +  CO     -->  CO    (MODE  I)  +  CO 

-11    -0.1293 
K  =  6.55x10        T  (Ref.  29) 


CO    (MODE  II)  +  H     — >  CO     +  H 

K  =  5.0xl0"14  (Ref.  29) 


CO    (MODE  II)  +  He  — >  CO    (MODE  I)  +  He 

K  =  5.895xl0"nT'°*1293  (Ref.  24 


CO     +  CO     —■ >      CO    (MODE  I)  +COz 

K=1.36xl0-17-808-17-8T"1/3  (Ref.  24) 
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CO    (MODE  I)  +  C02  -->  C02  (MODE  II)  +  C02 

,      .„        -11-0.1293    -1431. 5/T    ..  ,D    ,     0>n 

K=  6.55x10        T  e  vib  (Ref.  24) 
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APPENDIX  3:  HN~  +  C02  LASER  PROGRAM 
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